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Abstract. The dynamics of electrons in ac driven double quantum dots have 
been extensively analyzed by means of Floquet theory. In these systems, coherent 
destruction of tunneling has been shown to occur for certain ac field parameters. In 
the present work we analyze, by means of Floquet theory, the electron dynamics of 
a triple quantum dot in series attached to electric contacts, where the central dot 
position oscillates. In particular, we analyze the quasienergy spectrum of this ac 
driven nanoelectromechanical system, as a function of the intensity and frequency 
of the ac field and of external dc voltages. For strong driving fields, we derive, 
by means of perturbation theory, analytical expressions for the quasienergies of 
the driven oscillator system. From this analysis we discuss the conditions for 
coherent destruction of tunneling (CDT) to occur as a function of detuning and field 
parameters. For zero detuning, and from the invariance of the Floquet Hamiltonian 
under a generalized parity transformation, we find analytical expressions describing 
the symmetry properties of the Fourier components of the Floquet states under 
such transformation. By using these expressions, we show that in the vicinity of 
the CDT condition, the quasienergy spectrum exhibits exact crossings which can be 
characterized by the parity properties of the corresponding eigenvectors. 
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1. Introduction 

Triple quantum dot systems have been realized in the laboratory and their stability 
diagrams [T], charge [21 El S] and spin [5J [6] transport properties studied as a function 
of external parameters. These systems are important as components of quantum 
information circuits and molecular electronics. When one of the dots is set in motion by 
means of an external interaction, the system becomes a nanoelectromechanical system 
(NEMS). The study of these systems has become an important area of research, with 
potential applications as ultrasensitive detectors of mass [TJ, force [8], and displacement 
[9]. Recently, Azuma et al [10J have measured quantized tunneling current from an Au 
nanodot mounted on a cantilever, using a scanning tunneling microscope probe. Also, 
Koenig et al [11] reported the fabrication of a mechanical single electron transistor using 
ultrasound waves as the excitation mechanism. In all these works, a central role is played 
by an oscillating element acting as a shuttle carrying one electron from one contact to 
the other. The first proposal of such a device was made by Gorelik et al [T2] to be used 
as a single electron transistor. 

Ac driven quantum dot arrays have been extensively studied in recent years both 
experimentally p31 [HI [15] and theoretically [HI [17]. An important physical effect 
analyzed in these systems is the electron localization induced by an ac field, where 
for special parameter values of the driving field corresponding to roots of the Bessel 
functions [18], tunneling is suppressed. This effect, termed coherent destruction of 
tunneling, has been analyzed by means of Floquet theory and the corresponding 
quasienergy spectrum [19]. Close approaches of the Floquet quasienergies, as the system 
parameters are varied, produce large modifications of the tunneling rate and hence of the 
dynamics of the system. In particular, when two quasi-energies are close to degeneracy, 
the time-scale for tunneling between the states becomes extremely long, producing the 
phenomenon of CDT [2D] • 

In the present work, we extend the Floquet analysis [HH HTJ [221 [231 [21] to ac driven 
vibrating quantum dots arrays. In particular, we analyze the quasienergy spectrum of a 
Triple Dot Quantum Shuttle (TDQS), consisting of three quantum dots in series, with 
one level each, called Si, e c , and e r , respectively. The central dot oscillates with frequency 
u, and the system is in the presence of an ac- field characterized by frequency u ac and 
intensity V ac (figure [1]). Additionally, there is an infinite bias applied to the contacts 
of the triple quantum dot system, and there can be a voltage difference (detuning) 
Eb = si — e r between the left and right quantum dots. We refer to this setup as the 
dynamic case, to distinguish it from the static case where the central dot is fixed. For 
static three or more quantum dots in series, detailed studies (23, ES, EI] have been carried 
out looking for possible control of the electronic current features of the system from the 
analysis of crossings and anticrossings in the quasienergy spectrum. Here we study, both 
analytically, using first order perturbation theory, and numerically, the conditions for 
CDT in a TDQS, in terms of the quasienergies as a function of detuning, Eb, and the 
driving field parameters, in both static (a; = 0) and dynamic (u ^ 0) cases. For finite 
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detuning, in the static case, we obtain numerical results in agreement with a first order 
perturbative calculation in the weak tunneling regime, and find conditions where exact 
crossings of the quasienergies take place corresponding to CDT. In the dynamic case, 
we obtain different results depending on whether the current resonance as a function 
of detuning e& is calculated at even or odd multiples of the driving frequency u ac . For 
the former case, anticrossings are obtained at values of the driving field parameters 
corresponding approximately, as we will see below, to a root of a Bessel function. We will 
show how, in the case of weakly coupled quantum dots, these results can be reproduced 
with great accuracy considering interdot tunneling in first order perturbation theory. 
The resonances at odd multiples of u ac are shown to depend strongly on interdot tunnel 
coupling. For zero detuning, we note that the Floquet Hamiltonian is symmetric under a 
Generalized Parity Transformation (GPT), and the CDT phenomena can be explained 
by analyzing the quasienergy spectrum and the parity of the corresponding Fourier 
coefficients of the Floquet states. 

In section [2] we define the formalism used to describe our system with special 
emphasis in Floquet theory. We present our results in section [3l and finally the 
conclusions in section HJ 

2. Model 

We consider a linear array of three quantum dots where the central movable dot is 
flanked by two static dots at fixed positions ±x (figure [[]). The oscillation of the 
central dot affects the tunneling rates, which are now position dependent. The charging 
energy and interdot Coulomb interaction are assumed to be sufficiently large that only 
one transport electron can occupy the chain of three dots at any given time (Coulomb 
blockade regime). The TDQS is modelled by means of a tight-binding Hamiltonian 
H T qs = (H + H osc + Htun) where 

H = ei \l) (l\ + e c (x) |c) (c| + e r |r) (r| ; 

Hose — fojj d)d; 

H tun = T l (x)[\l)(c\ + \c)(l\]+ (1) 

T r (x) [|c> (r\ + \r) (c\). 

In the above equation, H Q represents the Hamiltonian of the three dots. Here the 
energy of the central dot e c depends on the position operator x given by x = Ax zp (cft+d), 
where Ax zp = {h/2muj) l l 2 is the zero-point uncertainty position of the oscillator. Thus 
e c (x) = [ei — (eb/2x )(x+x )], where it is assumed it undergoes a Stark shift proportional 
to its position due to the detuning across the device, Eb = (si — e r ); H osc is the 
Hamiltonian associated to the oscillation of the central quantum dot (frequency u>) 
where d\ d are creation and destruction boson operators respectively; H tun represents 
the tunneling between all three dots, with position dependent tunneling rates T r (x), and 
Ti(x), given by 

T,(x) = - Ve- a(xo+£) ; 
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Figure 1. Driven TDQS where oscillating tunneling barriers between the central 
quantum dot and its neighbors are considered. The system is driven by an ac-potential 
V(t) = V ac cos(uj ac t), and the left and right tunneling barriers are rigidly coupled to 
the leads with constant transition rates, T. The system is in an infinite bias setup, 
and there can be a voltage difference (detuning), Eh = £i — £ r , between left and right 
quantum dots. 



T r (x) = -V e - aix °- £ \ (2) 

where V is the tunneling amplitude, and a is the inverse of the tunneling length. The 
position operator x measures the displacement of the vibrational mode. The matrix 
elements of the tunneling couplings Tf 3 and T* J in the oscillator basis can be calculated 
analytically by using a procedure discussed in reference [28] which yields, 

Tf = - V [2"i\ j!] 1/2 e~ ax ° e (Q o /4) [-a ] i_i 

3 r 1 2k 

V IggJ i > i (3) 

^k\{k + i-j)\{j-k)V ~ h K) 

with T\ 3 = (—iy~ l Tp. In the above equation we have defined cto = (■s/2Ax zp a); the 
matrix elements for the case i < j are calculated by using the hermiticity of the tunneling 
couplings. Our model introduces a time- dependent ac-potential, V(t) = V ac cos(cu ac t), 
applied between the left and right quantum dots, and tunneling to and from the leads. 
Thus, the total Hamiltonian reads: H = H TQS + H ac + H ieads + H dot-leads, where 
H ac = (Vac/2) cos(u ac t) [\l) (l\ — |r) (r|] provides the effect of the ac-potential, which 
is introduced as an oscillation of opposite phase on the energy levels of the right and 
left dots flanking the central dot, and the contribution from the leads is given by 
Hieads = J2i J2k t e k c i ; c fcp where c 1 and c are creation and destruction operators for 
an electron with momentum k\ in lead /, respectively. Finally, Hdot-ieads stands for 
the coupling between the external dots and the corresponding leads, and is given by 
H dot -ieads = J2{i=L,R},k V ik c^ci + h.c, where is the dot-lead coupling matrix element, 
which leads to a coupling strength T; = 2n J2i % l^zfcl 2 — e ik)- 
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The equation of motion of the system using the standard master equation approach 
for the reduced density matrix (RDM) [29], generalized to include the environment 
of the oscillator, given by p^ s = [p^ s ]rdm + [p%]diss, where p^ s represents the matrix 
element with the dot states q, and s (q,s — 0,/,c, r), with the oscillator states i, and 
j = 0,1,2,...). The term [p^ s ]rdm incorporates transitions between the leads and 
the outer dots, 

( T d^ q Pd d - T q ^ d p%) ;q = s 
Ed*, T i^d + Ed^ s T s^d) P%\ q^s. 

Here we have also considered the transitions r r _>.o — To-*/ = T, and the term 
[Pq j s ]diss = -7[(«+j)/2)]p^+7[(« + l)(j + l)] 1/2 Pgt 1 ' j+1 , accounts for the dissipative effects 
of the oscillator's environment [30] in the low temperature limit; here 7 stands for the 
classical damping rate of the oscillator. The calculations for the current are performed 
by numerical integration of the RDM and averaging over the resulting electronic current 
measured at the right quantum dot, i.e. I = eT ^ [p"rW 

2.1. Floquet theory for oscillating triple quantum dots 

Since our system is described by a time periodic Hamiltonian, we can use Floquet 
theory [2H [22] in order to analyze the quasienergy spectrum of the system. We 
obtain the quasienergies by numerical diagonalization of the Floquet Hamiltonian 
Ti = [H' — ihd/dt], (with H' = Htqs + H ac ), represented in the composite basis 
(£) O T of quantum dots (]/), |c), \r)), oscillator = 0, 1, 2, ...N) and the periodic 
functions of t with period r = (27r/c<; ac ), respectively. The matrix elements in this 
composite representation are given by 

U% n ' n = ((n'\H%\n)} = - / dte- in ' u *° t H%e inu '»> t , (5) 

T Jo 

where we have used the definition of the inner product in space T, involving a time 
average. If N t is the number of photon sidebands, then n, n' = 0, ±1, ±2, ±3, ±N t . 
This yields, 

H]f n n = (e t + jhu + nhuj ac )5 n i n 5i 



Vac, 

4 



H — (8n',n+l + fin',n-l)fiiji 



UT' n = - Ax zp -^{^/j + l5 iij+1 + sfj 8^)8^ 

ZXq 

+ (jhu + nflWa^SijSn'n; 
W r j ; n ' n = (e r +jhhJ + nhlOa^n'Jij 

—{3ri,n+l + $ri ,n-l)&ij] 

njij,n'n _ njij,n'n _ rpij/£\c 

tt lc - n d - i l [x)o n > n , 
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n lr - H rl - u; 

nT' n = ni j f n = T^(x)5 nln . (6) 

The Floquet quasienergy spectrum of the system can be obtained by numerical 
diagonalization of the matrix whose elements are given above. The eigenvectors of T-L in 
the composite basis P(^)00T can be written as 

\Mq;j,t)) = J2 c n q ^- l7mJ act W,j), (7) 

n,q,j 

where the kth Floquet state, is represented in the composite basis \q;j) = \q)(£)\j), 
where the indexes q and j correspond to the dot and oscillator states, respectively 
Thus, the Floquet eigenstates are given in effect by the Fourier coefficients C* •. As we 
shall later see, the parity of the Floquet states determines the symmetry properties of 
these Fourier coefficients. 



2.2. Perturbative approach 

In this section we describe an alternative method for exploring the quasienergy spectrum 
for both static and dynamical systems. In both cases, we consider that the tunneling 
Htun can be treated as a perturbation, i.e., T <C V ac [16j [31], allowing the application 
of standard Rayleigh-Schrodinger perturbation theory. We present here the details of 
the perturbative treatment for the dynamic case and the static case is given in the 
Appendix. 

We consider for simplicity the case N = 1 in the oscillator basis (i, j = 0, 1). We 
start from the Floquet equation corresponding to the Hamiltonian 

[H , -ihd/dt)\4>{t)) = e i \<l>(t)), (8) 

with H' = (Htqs + H ac ). The unperturbed Floquet Hamiltonian given by H u = 
(Ho + H ac + H osc — ihd/dt) is approximately diagonal in the composite basis V(£)0 
of the quantum dots and the oscillator (whenever (Ax zp Eb/2x Q ) <C Hlj), with matrix 
elements U% qs = (i\H u ,qs\j) , given by, 



n/ij 




(si + jhw) + ^cosu ac t 




Oij] 


H ij 

' *~u,cc 




d~ 

ei - £ b /2+jhu - ik- 
on, 






H ij 




(e r + jhu) - ^cosu ac t - ih^ 




n/ij 


- ^u,d - °; 






n/ij 


_ qjij — n- 
— n u,rl ~ u > 






H ij 


= = 







with i,j = 0,1. For the particular case of a symmetrical TDQS the unperturbed 
Hamiltonian H u is exactly diagonal in the composite basis T> (^) O, and no 
approximations are needed. 
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Following the same procedure as in the static triple dot case discussed in the 
Appendix, we obtain the following eigenvectors, 



Ol) 


= e~ 


-i(e b /2~t Q i)t/h e ~i( 


sino; ac t|g j^y 


I) 


= e~ 


-i(e b /2-e u )t/h e -ii 


sinw ac i| -y j^y 


Oc) 


= e~ 


-fcoct/»| 0c ) ; 




1c) 


= e~ 


-fclct/»| lc ^ ; 




Or) 


= e~ 


-i(e b /2-e 0r )t/h -it 


•sinw ac t|Q r \. 


I) 


= e~ 


-i(e b /2~ei r )t/h g — il 


Jsinaj ac t|^ ^\ 



(10) 

where £ = V ac /2u ac . 

The corresponding Floquet eigenenergies are given by, 

e ol = Si + n (0) hu ac ; 

e u = Ei + ku + n {1) hu ac ; 

eoc =ei- (e 6 /2) + n {3) frw ac ; 

e lc = Si - (e b /2) + hu + n {4:) ftw ac ; 

e 0r = E r + n i5) hu ac ; 

€\r — e r + hljJ + n^huj ac , (11) 

where n (fc) = 0, ±1, ±2, with k = 1, 6. 

We now proceed to calculate the matrix elements of the perturbation Htun m 
the basis that is, we compute P % q { = ((<j>f <a \Ht un \<f>j tS )), where i,j 

q,s = l,c,r. We can verify that P\f 



qs 
pij 



pi'j _ p l J _ p l J 
rr lr rl 



0, 1, and 
0, while the other non- 



zero matrix elements can be calculated by using the identity given by equation flA.llj) 



and by considering two important conditions: (i) that an alignment of the ac-sidebands 



occurs i.e. €j C = e^, and e JC = e 3 >, whenever £ b = 2s'hco ac , (as in the static case), and 



that (ii) hjj = n'hu ac (where s' and n' are positive integers); the 6x6 matrix P = (P* J S ) 
is then given by 



{Pa) (Pic 

(Pel) (Pec 
(Prl) (Pre 



(Plr) 
(Per) 
(Prr) 



(12) 



where the 2x2 submatrices are 

{Pic) - 



7f J. 



(0 r^j^AO 



(Pel) 
(Per) 



t™J-A0 T^j^Ai) 

T l J n '- s >(C) Tl l J_ s i(£) 

T™JAt) T r 01 J s ,_„,(0 " 

T™j n , + Ai) T r ll JA0 



(13) 
(14) 
(15) 



Quasienergy spectrum and tunneling current in ac-driven triple quantum dot shuttles 8 



(Pr, 



T™J s ,(0 T r m J n , +s ,(0 



(Pi) = (Pec) = (Prr) = {Plr) = (Pi 








(16) 



(17) 



The quasienergies are computed by a numerical diagonalization of P (equation (fl2l)); in 
the case of the symmetrical (e^ = 0) TDQS, with u ac = 0.2, and u — 1, we get s' — 
and n' = 5. For the asymmetrical case (e& 7^ 0), and the particular case of Eb = 0.4, we 
have s' = 1 in equation f fT2l . 

For the static case (H osc = 0) we obtain the following analytical expressions for the 
quasienergies (the details are discussed in the Appendix), 

"£b _ 
.2 

(— ei) (mod fkjj a 



V2T J S '(0] (mod huj, 



(mod TujJn 



:ic c 



where T = — Vexp(— axo). This result is in agreement with a general perturbative 
treatment of ac-driven tight-binding chains of quantum wells or quantum dots [251 [26J 
127] . Furthermore, it can be seen that an exact crossing between all three quasienergies 
is obtained whenever £5 is an even multiple of hu ac and it corresponds to the roots of the 
corresponding Bessel function J s > (See for example figure H] (b), where, for the chosen 
parameters, s' = 1). 



3. Results 



In what follows, we shall explore the main features of the current at finite detuning, 
and the quasienergy spectrum for the TDQS driven by an ac gate voltage in the weak 
interdot coupling regime, with special emphasis on the CDT phenomenon. All of the 
calculations are worked out in units where h = 2m = e = 1. In figure [2] (a) we present 
the electronic current, J, as a function of detuning Eb, where a series of resonances 
at even and odd multiples of the ac-frequency u ac can be seen. In a previous work 
[32J, we analyzed numerically the tunneling current through a Triple Quantum Dot 
(TQD) oscillating device and found that the positions of these peaks obeyed sum rules 
indicating the number of photons emitted or absorbed in the process (v) and the final 
oscillator state (n) ((v,n) for abbreviation). In figure [2] (a) we present results for the 
current as a function of detuning e&, for two values of the intensity of the ac-field V ac . 
For V ac = 1 (black dashed line), we note that there are current resonances at values 
of Eb which are integer multiples of u ac . Stronger peaks are observed at even multiples 
of u ac corresponding to the alignment of ac-sidebands developing from eigenstates q, e c 
and e r . The weaker peaks at odd multiples of u ac , correspond to the alignment of levels 
€1 and e r . These are due to the overlap of left and right quantum dots, as discussed in 
Ref. 132]. 
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Figure 2. I versus Sb for a triple dot system for the ac parameters u! ac — 0.2, V ac = 1 
(black dashed line) , V ac = 1.533 (red solid line), corresponding to first root of Ji(£), 
and typical values of the parameters {V; 7; xo] T; a}={0.5; 0.01; 5.0; 0.05; 0.4}, for 
both the (a) dynamic (w = 1 = 5uj ac ), and (b) static (ui = 0) cases. In (a) we also 
include for comparison the undriven current (V ac = 0, blue dotted line). Note the 
reduction of the resonance peak at e& = 0.4 (red solid line) in both cases, and the 
appearance of resonances at odd multiples of uj ac in (a) . These are due to the overlap 
between wave functions of left and right quantum dots as discussed in Ref . [32! . 

Let us recall under what conditions the electronic current exhibits the CDT 
phenomenon. The current peak at e b = 0.4 (and u ac = 0.2), corresponds to v — 1 
emitted photons and n = is the final oscillator state (see Ref. [32J). Interestingly, 
if we calculate the current for £ = 3.832 {V ac = 1.533, red solid line in figure [2] (a)) 
corresponding to the first root of with all other parameters fixed, we see that the 

peak at Eb = 0.4 decreases significantly. This is the phenomenon of coherent destruction 
of tunneling in our triple dot quantum shuttle device. Additionally, we note that at 
Eb = 2.0, the opposite situation occurs, where the resonance for V ac = 1 (dashed black 
line) is small compared to the case V ac = 1.533 (red solid line). The reason is that this 
peak corresponds to v = 0, n = 1 and V ac = 1 is very near the value corresponding to 
the first zero of Jo(0 (Vac = 0.96192, corresponding to £ = 2.4048). Included in this 
figure for comparison purposes is the current for the undriven case (blue dotted line), 
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Figure 3. / versus Eb for a triple dot system in the dynamic case for a smaller 
tunneling than in the previous case of figure [5] (a), by choosing a — 0.8 with ac 
parameters w ac = 0.2, V ac = 0.5 (black dashed line) and V ac — 0.96192 (red solid line), 
corresponding to first root of Jo(£) (£, = 2.4048), and values of the parameters {V; 7; 
x ; r}={0.5; 0.01; 5.0; 0.05}. Also shown is the current for the undriven case (blue 
dotted line). We can see that the undriven current exhibits two well defined resonances 
at Eb = 0, and 2.0. These resonances are effectively turned off by adequately tuning 
the field intensity V ac to a CDT condition. 

where resonances appear at multiples of the oscillator frequency u = 1. In figure |2] (b), 
we present the results for the static case (co = 0). Here we see that current resonances 
appear only at values of e b which are even multiples of u ac . These resonances correspond 
to v — 1, 2, 3, ... emitted photons involving the final oscillator state n = as expected. 
Note that the peak at e b = 0.4 practically disappears at the CDT condition (red solid 
line). Thus, for this particular choice of parameters, the CDT phenomenon is more 
evident in the static case than in the dynamic case, where additional tunneling channels 
(besides v = 1) contribute to the current. This will be corroborated further when 
we analyze below the properties of the quasienergy spectrum. Interestingly, we also 
note the current peaks at e b = 1.6 and e b = 2.0 for V ac = 1 (black dashed line) have 
also decreased. According to the sum rules, these peaks correspond to v = 4 and v = 5 
emitted photons, respectively, i.e. interdot tunneling is renormalized by Bessel functions 
J 4 (0 and J 5 (0 which, for V ac = 1,£ = 2.5, are small (< lO" 2 ). 

Next, in figure [3] we plot the electronic current I as a function of detuning e b in a 
TQD oscillating system (i.e. dynamic case with oj = 1) for smaller tunneling couplings 
Ti and T r (i.e. smaller than in the case of figure 2 (a)), by choosing a value of a = 0.8. 
The undriven current (blue dotted line) exhibits two strong shuttling resonances at 
values of s b = 0, and 2.0, which correspond to alignment of the energy levels of the 
vibrational subbands of the three dots. When the ac voltage is applied (black dashed 
line), an additional energy level splitting occurs which gives rise to the new resonances 
at even multiples of u ac . By an adequate manipulation of the driving field intensity, 
we can induce the CDT condition (red solid line) in order to turn off the shuttling 
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Figure 4. Static case: (a) Quasienergy versus e& for ac parameters Lo ac — 0.2, and 
V ac — 1 (black hollow dot line) and V ac = 1.533 (red solid dot line), and typical values 
of the parameters {V; x ; T; a}={0.5; 5.0; 0.05; 0.4}. Note the crossings of the red 
line at e = 0.4, and of the black line at e = 1.6, and 2.0, both indicative of CDT. (b) 
Quasienergy versus £ = (V ac /2u; ac ) for e = 0.4, and same parameters as in (a) (red 
hollow dot line). Results from a first order perturbation calculation given by equation 
(jlH (blue solid line), are also shown. 



resonance peaks corresponding to E\> = 0, and 2.0. Also note there are no resonances at 
odd multiples of the driving frequency. The reduction in tunneling effectively cancels 
the overlap between functions of left and right quantum dots. 

3.1. Quasienergy spectrum analysis 

Next, we carry out a numerical calculation of the quasienergies as a function of detuning 
6b and driving field parameter £. Following Fromherz [23J and Je et al [34], we use 
N t = 20 Fourier components in order to get accurate results for the range of values of 
the driving field intensity V ac . 

In the static case, when we examine the quasienergies as a function of detuning 
£b (figure H] (a)) we see the three quasienergies ei,e c ,e r , along with their replicas at 
multiples of u ac (V ac = 1-0, black hollow dot line), in agreement with the Floquet 
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Figure 5. Dynamic case: (a) Quasienergies as a function of detuning e b for V ac = 1 
(black solid dot line) and V ac = 1-533 (red solid dot line), corresponding to CDT 
conditions for Sb — 0.4, with io ac — 0.2 and typical values of the parameters {V; 
xo] T; a}={0.5; 5.0; 0.05; 0.4} with N t = 20. (b) Quasienergies as a function of 
£ = (14 c /2w ac ) for £{, = 0.4 with ui ac — 0.2 and the rest of the parameters the same 
as in (a) (black solid dot line). The result of a first order perturbative calculation 
obtained from equation (fl"2j) (red hollow square line) is included for comparison. 



theorem. Furthermore, we see anticrossings involving all three levels at the values of 
e b where current resonances appear. For V ac = 1.533 (red solid dot line) corresponding 
to the first zero of Ji(£) (CDT condition), an exact crossing develops at e b = 0.4, in 
agreement with the disappearance of the current resonance at this detuning (figure [2] 
(b)). Two more CDT conditions occur at e b = 1.6 and e b = 2.0, indicated by the 
crossing of the black hollow dot line at these values, as discussed before. Also for the 
static case, in figure @] (b) we show the quasienergies as a function of £ for the case 
e b = 0.4. Note the exact crossings (CDT) at £ ~ 3.832 and 7, corresponding to the first 
two roots of the Bessel function Ji(£), and to the quenching of the current resonance 
at this detuning in figure 2 (b). We see that the numerical results agree quite well with 
the perturbation calculation (equation (TlBl). using Ji(£). 

For the dynamic case, in figure |5] (a) we show results for the quasienergies as a 
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Figure 6. Quasienergies as a function of £ = (V QC /2cL> ac ) for e& = with u> a c = 0.2 
and typical values of the parameters {V; x^; r}={0.5;5.0;0.05}, with N t — 20 (black 
solid dot line), for the case a = 0.8 of figured] Note the excellent agreement obtained 
as compared with the result of a first order perturbative calculation obtained from 
equation (TT2"|) (red hollow square line), where both plots almost completely overlap. 



function of detuning e b , for the same two values of the ac-field intensity. For V ac = 1 
(black dot line) we see a pattern of crossings (anticrossings) where Sf, is an odd (even) 
multiple of u ac where current resonances appear. Interestingly, for V ac = 1.533 (red 
line in figure [5] (a)), corresponding to the CDT condition, we note a drastic reduction 
in the anticrossing separation (compared with the previous case V ac = 1) at £& = 0.4 
(the region highlighted by the blue square in figure [5] (a)), which signals the appearance 
of CDT. Thus, even for the case of finite detuning, where no symmetry properties can 
be invoked for the eigenstates, the presence of CDT can still be identified with these 
near crossings in the quasienergy spectrum. This is further corroborated by plotting the 
quasienergies as a function of £ (figure [5] (b)), for fixed Eb = 0.4, where again we obtain 
these near crossings at the roots of the Bessel function Ji(£). Also plotted in this figure 
are the results of the perturbation calculation, equation ( TT2l . using Ji(£) (red hollow 
squares). Note that the perturbation calculation also reproduces the near crossing. In 
figure [6] we show the quasienergies as a function of £ (black solid dot line), for the case 
of a smaller value of tunneling (a = 0.8) compared with the case discussed in figure 
El for Eb = (symmetrical system), and u ac = 0.2. In this case the CDT phenomenon 
occurs at a particular value of V ac = 0.96192, with a ratio £ = 2.4048 that corresponds 
to the first zero of Jo(0- The perturbation calculation using equation f fT2l (red hollow 
square line) gives an excellent agreement with the numerical results. 

For completeness, let us explore more features of the CDT phenomenon in TQD 
oscillating symmetrical systems (Eb = 0), by considering the case presented in figure [5] 
(a). We perform a Floquet analysis of the CDT by analyzing the quasienergy spectrum 
of the system. For our calculations, we choose N = 1, i.e. two oscillator states 
n = 0, 1 for each quantum dot. In figure [7] we plot the quasienergies as a function 
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Figure 7. The inset shows the quasienergies in the first Brillouin zone (delimited with 
blue solid lines to help the eye) as a function of £, for Sb = 0, with fixed oj ac = 0.2, 
N = 1 oscillator states, N t = 20 and the rest of the parameters the same as in 
figure [5] (a). The main graph is an amplification of the quasienergies in the vicinity of 
£ = 2.4048 (indicated by a blue arrow) which corresponds to the zero of Jo(£) where 
CDT occurs. The regions A, and B (see text) exhibit two crossings of the quasienergies. 
The quasienergies are classified according to parity i.e. even states (red solid dot line) 
and odd states (black hollow dot line). 



of the ratio £. The inset of figure [7] shows the quasienergies in the first Brillouin zone 
(—cj ac /2 < e a < ui ac /2) as a function of £, where we can appreciate a complex pattern 
of crossings and anticrossings of the quasienergies in the vicinity of £ = 2.4048; an 
amplification of the region where this behavior occurs is illustrated in the main graph. 
Here we can identify regions A and B, where two crossings of the quasienergies take 
place. The corresponding eigenstates exhibit opposite parities, where the even states are 
denoted by red solid dots and the odd states with hollow black dots. The above results 
are supported by detailed numerical studies, and by using an interesting symmetry 
property of the Floquet Hamiltonian. At zero detuning, the Floquet Hamiltonian is 
invariant under a generalized parity transformation Sqp '■ (x,t) —> (—x,t + r/2) |24| 133]. 
Hence, the Floquet states can be classified according to parity. From equation (J7|), we 
can write the Fourier coefficients, ... as 

C Lj = - [ {q;m)j™ aot dt, (19) 

and, after performing the Sgp transformation, it can be shown that the above coefficients 
obey the following relations, 

c^ j = ±(-ir(-iycl_ q>j , (20) 

where the upper (lower) sign stands for even (odd) states, k stands for the k-th Floquet 
state, and the indexes are defined as before: n, q, and j, associated to the temporal, dot, 
and oscillator basis, respectively. In our notation, the — q index in the right hand side of 
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the above equation indicates a spatial site inversion of the dots, i.e. I -H- r. Using these 
properties, we have verified that at A and B (figure [7]) the corresponding eigenstates 
exhibit opposite parities. 

4. Conclusions 

In summary, we have carried out a Floquet analysis in an oscillating ac-driven triple 
quantum dot system in order to investigate the quasienergy spectrum as a function 
of the field parameters and the external voltage. The CDT condition is discussed in 
terms of the properties of the quasienergy spectrum as a function of detuning and ac 
parameters. For that purpose, in the strong driving regime, we extend a perturbative 
analysis for the static case to the dynamic configuration. For finite detuning, we find 
conditions under which, CDT is realized in both static and dynamic cases. In the 
latter smaller reduction of the current resonance peak is obtained as compared 

with the static case, due to additional channels contributing to the tunneling current. 
We have also found in the dynamic case that the ac drive, for certain field parameters, 
induces a CDT condition that destroys the tunneling current peaks (tunneling channels) 
for the undriven oscillator. Therefore, ac fields are able to modify and even block the 
tunneling current otherwise flowing through an oscillating triple quantum dot. From 
the analysis of current-voltage characteristics, we obtain current resonances at voltages 
which are integer multiples of the driving frequency u ac - For the case of even multiples 
of u ac , in the static case, the numerical results for the quasienergies can be reproduced 
by analytical perturbative calculations. Exact crossings of the quasienergy spectrum 
are obtained, characterizing the CDT condition, while near crossings are obtained in 
the dynamic case. For the case of odd multiples of u ac , the numerical result cannot 
be reproduced by first order perturbative calculations. With respect to this case, it 
is important to stress that interdot tunneling is coherent in the present model and, 
therefore, tunneling is included to all orders in our numerical calculation. However, this 
is not the case for the perturbative model presented here, where we have only considered 
first order perturbation theory. For zero detuning, by exploiting the invariance of the 
Hamiltonian of the system under a generalized parity transformation, we have found 
that the Fourier components of the Floquet states satisfy analytical expressions which 
take into account the parity of the oscillator states. These allow us to classify the 
quasienergy spectrum of the oscillating driven TQD according to parity. By using these 
expressions, we show that, at zero detuning, in the vicinity of the CDT condition, the 
quasienergies exhibit crossings, where the eigenvectors are characterized by opposite 
parities. Our results are a contribution in ongoing efforts to control electron transport 
in few electron nanoelectromechanical quantum dot arrangements. They should be 
relevant, for example, in Molecular Electronics, or in the area of quantum information 
circuits where these systems have been considered as spin entanglers [35] or quantum 
switches |36|. 
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Appendix: Perturbative approach for the static case 

In this section we seek analytical expressions for the quasienergies using a perturbative 
approach for the static triple dot. We start from the Floquet equation corresponding to 
the Hamiltonian, 

[H" - iW/dt) \<pi(t)) = €i |&(t)>, (A.l) 

with H" = (H + H tun + H ac ). 

In this approach the tunneling H tU n is treated as a perturbation, i.e. T <C V ac 
[16] . using Rayleigh-Schrodinger perturbation theory. This procedure deals with the 
diagonalization of the perturbation matrix in the base of the unperturbed Floquet 
Hamiltonian % u given by H u = (H + H ac — ihd/dt). 

We proceed to calculate the eigenvectors and eigenvalues of H u starting from the 
equation, 

= e?|#(t)>. (A.2) 

The matrix elements of the unperturbed Hamiltonian in the dot basis T> i.e. 
7-Lu, qs = (q\Hu\s) (q, s = l,c, r) are given by, 

Hu,u = — + — cosu ac t - %h—\ 
2 2 at 



at 





£b 


Vac 


u,rr 


2 


2 


u,lc 


— Wu,cl 


= 0; 


u,lr 




= 0; 


u,cr 


T~Lu,rc 


= 0, 



cosu ac t - ih—\ 

at 



(A.3) 

The eigenvectors can be explicitly calculated from the matrix representation of equation 
f)A.2p . and to exemplify the procedure we consider the particular case for <fif(t), which 
leads us to the equation, 

(| - e?) + (Kc/2) cosu ac t\ tf(t) = ihd<ft(t)/dt, (A.4) 

which has the solution 

0«(t) = e -*(^/2- e r)*/» e -««T««act j ( A-5 ) 
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with £ = (y ac /2hu ac ). The corresponding eigenvector reads, 

Similarly, we also obtain, 



and 



e —i(—£l,/2-e%)t/h e +«? sinuiact 



r) 



(A.7) 



\c) . (A.8) 

The eigenvalues are calculated using the periodicity of the Floquet states i.e. (<fif(t) = 
(j)f(t + 2ir/u ac )) and are given by ef = (eb/2 + m^huiac), = + m^Huj ac ) and 

e « = (0 + mWhhj ac ), where m< fe) = 0, ±1, ±2, for fc = 1, 2, 3. 

We now proceed to calculate the matrix elements of the perturbation H tun as 
follows, 

1 



qs 



H 



tan 



dt 



H 



tun 



(A.9) 



where q,s = I, c, r. From the above equation it can be readily verified that Pu = P cc = 
P rr = P\ r = P r i = 0, and the remaining non-zero matrix elements can be calculated 
by considering that an alignment of different Floquet quasienergy ladders occurs i.e. 
e" = e", and e" = ef, whenever the condition e h = 2s'Hcu ac (s' = 0, 1,2, ...), is satisfied. 
For the particular case of Pi c , 



Pic 



Htun 

dt 



(A.10) 



o 



H 



tun 



T 1 - 

T 



Pic 

by using the identity 

± it; sincoact 



^j. e +ii sinu) ac t ^Eo/2-e^t/h ^ie™ t/h 



we obtain 



Pic = T £ J B (0- / 



^ ^iulacS't ^iuj ac nt 



(A.11) 



(A.12) 



which leads us to p c = T J- s '(£) = — T J s '(0- Similarly, one can show that P C £ = P c . 
P rc = —TJgi(£), thus the matrix elements of the perturbation are given by, 



P = -TJ s ,{0 



1 

1 1 
1 



(A.13) 



and the diagonalization of equation (IA.13|) leads us to the set of eigenvalues 
{0, ±y/2T J s '(£)}, where the result given by equation f[T%j) follows. 
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